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Abstract 

We study two-dimensional turbulence in a doubly periodic domain driven by a 
monoscale-like forcing and damped by various dissipation mechanisms of the form 
z/^(— A)^. By "monoscale-like" we mean that the forcing is applied over a finite range 
of wavenumbers k^in < k < fcmaxj and that the ratio of enstrophy injection 77 > to 
energy injection e > is bounded by k'^^^e <rj< k'^^^e. Such a forcing is frequently 
considered in theoretical and numerical studies of two-dimensional turbulence. It is 
shown that for > the asymptotic behaviour satisfies 
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where ||u|| and \\u\\i are the energy and enstrophy, respectively. If the condition of 
monoscale-like forcing holds only in a time-mean sense, then the inequality holds in 
the time mean. It is also shown that for Navier-Stokes turbulence = 1), the time- 
mean enstrophy dissipation rate is bounded from above by 2i'ik'^^-^. These results 
place strong constraints on the spectral distribution of energy and enstrophy and of 
their dissipation, and thereby on the existence of energy and enstrophy cascades, in 



Preprint submitted to Elsevier Science 



8 February 2008 



such systems. In particular, the classical dual cascade picture is shown to be invalid 
for forced two-dimensional Navier-Stokes turbulence (/i = 1) when it is forced in this 
manner. Inchision of Ekman drag {fi = 0) along with molecular viscosity permits a 
dual cascade, but is incompatible with the log-modified —3 power law for the energy 
spectrum in the enstrophy-cascading inertial range. In order to achieve the latter, it 
is necessary to invoke an inverse viscosity < 0). These constraints on permissible 
power laws apply for any spectrally localized forcing, not just for monoscale-like 
forcing. 
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1 Introduction 



It has long been recognized (see for example Fj0rtoft [9]) that the simultaneous 
existence of two quadratic inviscid invariants, energy and enstrophy, drasti- 
cally changes the picture of two-dimensional (2D) Navier-Stokes turbulence 
in comparison with that of its 3D counterpart. In the 1960's it was suggested 
by Kraichnan [11], Leith [14], and Batchelor [2] (hereafter referred to as KLB) 
that these two constraints give rise to the realization of two distinct inertial 
ranges in wavenumber space for forced 2D turbulence. More precisely, the KLB 
theory predicts that for a 2D Navier-Stokes fluid driven by a spectrally lo- 
calized forcing, energy is transferred to smaller wavenumber s while enstrophy 
is transferred to larger wavenumbers. For a system of finite size the energy is 
then predicted to cascade to the long- wavelength end of the spectrum up to the 
largest scale available (the system's hnear length scale), while the enstrophy is 
predicted to cascade to the other end of the spectrum, down to a dissipation 
length scale. This prediction is well confirmed numerically in transient evo- 
lution from spectrally localized initial conditions. An immediate corollary of 
the dual cascade and the hypothesis of a scaling symmetry is that the energy 
spectrum scales as k~^/^ in the energy-cascading range and as (with a 
possible logarithmic correction) in the enstrophy-cascading range, where k is 
the wavenumber. The latter scaling is consistent with the hypothesis that the 
dissipation of enstrophy is confined to the small scales; however, the former 
scaling is inconsistent with the hypothesis that the dissipation of energy is 
confined to the large scales. Hence, these scaling laws are incompatible with 
a persistent dual cascade in finite systems. This suggests that there is a prob- 
lem with the KLB theory applied to finite systems. (We refer here only to 
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the forced-dissipative equilibrium behaviour, not to transient evolution from 
spectrally localized initial conditions.) 

For the case of a time-independent force, Constantin, Foias and Temam [7,8] 
have proven the existence of a compact global attractor for forced 2D Navier- 
Stokes turbulence in a bounded domain. Constantin, Foias and Manley [6] 
have furthermore shown that for the special case of a doubly periodic domain 
and forcing of a single length scale, the KLB scaling laws cannot be achieved 
on the global attractor. Although this strong result appears to prove the un- 
realizability of the KLB theory for forced 2D Navier-Stokes turbulence in a 
finite domain, as argued heuristically in the previous paragraph, the question 
arises whether it is an artifact of the special choice of a constant monoscale 
forcing. 

In this paper we extend the results of Constantin, Foias and Manley [6] in 
several directions. First, for their choice of a constant monoscale forcing we 
derive an asymptotic bound on the enstrophy in terms of the energy and show 
that the time-mean enstrophy dissipation occurs around the forcing scale, 
precluding the existence of an enstrophy-cascading inertial range (whatever 
the energy spectrum). This result applies for any dissipation operator of the 
form A)'* with // > 0, not just molecular viscosity. Second, we show that 
the result applies for any "monoscale-like" forcing, by which we mean forcing 
over a range of wavenumbers /cmin < k < kmax where the ratio of enstrophy 
injection 77 > to energy injection £ > is bounded by k^i^^s < 77 < k^^^e 
for all possible velocity fields u. This is a classical (though not exclusive) 
scenario for the KLB theory (e.g. Kraichnan [11], p. 1421b; Pouquet et al. 
[21], p. 314; and Lesieur [15], p. 291), and is furthermore a common set-up in 
numerical simulations of forced 2D turbulence (e.g. Lilly [16]; Basdevant et 
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al. [1]; Shepherd [22]). Qlf the condition of monoscale-hke forcing holds only 
in a time-mean sense, then the time-mean results still go through and the 
asymptotic bound is replaced by a bound involving the time-mean energy and 
enstrophy. 

The question then arises whether the KLB theory and the dual cascade picture 
can be recovered, for monoscale-like forcing, with other dissipation mecha- 
nisms. It is shown that the introduction of Ekman drag (yU = 0) — a commonly 
used numerical device and one with some physical justification for geophys- 
ical applications [20] — in addition to regular viscosity permits the existence 
of an energy dissipation range at large scales and an enstrophy dissipation 
range at small scales, thus allowing a dual cascade. This is consistent with 
numerical results, e.g. Boffetta et al. [3], which demonstrate an inverse en- 
ergy cascade and a —5/3 power law under these conditions. However, the 
log-modified —3 power law of Kraichnan [12] for the energy spectrum in the 
enstrophy-cascading range is shown to be unachievable; the spectrum must be 
algebraically steeper than —3. This result is consistent with numerical sim- 
ulations of forced 2D turbulence under these conditions, which find steeper 
spectra, (e.g. Maltrud and Vallis [18] find the spectral slope of the enstrophy- 
cascading range to be between —3 and —4: the majority of their simulations 
yields values between —3.3 and —3.6. They also recover the scaling k~^/^ in 

^ Constantin, Foias and Manley [6] showed that the KLB scaling is potentially 
realizable (not to say it actually occurs) if the forcing simultaneously injects en- 
ergy at higher wavenumbers and removes energy at lower (though possibly nearby) 
wavenumbers. In this case the ratio of enstrophy to energy injection can be much 
larger than the square of the characteristic forcing wavenumbers and could actually 
lie in the dissipation range. Such a forcing is not monoscale-like by our definition. 
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the energy-cascading range.) 



It is finally shown that allowing an inverse viscosity (// < 0) together with 
regular viscosity does permit KLB scaling. It is notable that the numerical 
simulations of Borue [4], and the more recent high-resolution numerical sim- 
ulations of Lindborg and Alvelius [17], which both claim to exhibit a —3 
power-law enstrophy-cascading range, employ an inverse viscosity. The rather 
surprising implication of our results is that verifying the KLB theory in the 
case of monoscale-like forcing does not only depend on achieving sufficiently 
high Reynolds numbers (as is commonly believed), but also depends on the 
nature of the rather ad hoc dissipation operator employed at the large scales. 

The alternative, considered by Constantin ct al. [6], is to abandon the notion 
of monoscale-like forcing. This interesting possibility is not considered here. 

The remainder of this paper is organized as follows. Section 2 describes the 
2D Navicr-Stokcs equations and their mathematical setting, together with 
some basic inequalities. The concept of monoscale-like forcing is discussed 
further, and some examples given. Section 3 reviews some classical ideas and 
arguments in 2D turbulence. After those preliminaries, section 4 reports an 
asymptotic analysis, the results of which include a dynamical constraint and 
an upper bound on the energy dissipation rate. Section 5 extends the results 
of section 4 in an attempt to bound the time-mean enstrophy dissipation rate 
and explore possible scaling laws for the energy spectrum. These results in fact 
apply for any spectrally localized forcing, not just for monoscale-like forcing. 
The paper ends with some concluding remarks in the final section. 
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2 Governing equations and basic inequalities 

We consider 2D incompressible fluid motion confined within a doubly periodic 
rectangular domain Q. The fluid is assumed to be driven by a monoscale-like 
forcing (to be discussed further below) and damped by a variety of possible 
dissipation mechanisms including Ekman drag, hypoviscosity, molecular vis- 
cosity, hyperviscosity, and inverse viscosity. The 2D Navier-Stokes equations 
which govern the fluid motion are written in abstract form in a function space 
H as 

(All 

— + B{u,u) + Y.^,A'^^^f, (1) 

u{t = 0) =Uo, 

where > is a generalized viscosity coefficient, A = —A, and / is the 
forcing. The number // will be called the degree of viscosity. When n — 1 we 
have the usual molecular viscosity, while = corresponds to Ekman drag. 
The cases in which /i > 1 and < /x < 1 correspond to hyperviscosity and 
hypoviscosity, respectively. We also consider inverse viscosity, i.e., negative 
values of /i. The summation is taken over all dissipation channels involved. 
A detailed description of the functional analysis setting for (1) is given in 
Constantin and Foias [5] or Temam [25,26]. We recall that H is the L^-space of 
periodic, non-divergent functions with vanishing average in fl. B{u, u) = P{{u- 
V)ii) where P is the orthogonal projection in onto H. We denote by if" the 
domain of definition of A"/^ for real a. The (degenerate) positive eigenvalues 
of A are denoted by with the index k being the wavenumber, and the 
eigenspace corresponding to is denoted by H{Xk). We will occasionally refer 
to A^^^^ (or simply k~^) as the length scale associated with the wavenumber k. 
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The scalar product and the norm in if" are given respectively by 



[u 



,v)a = Ju-A''vdx, (2) 

\\uL-(u,uyj\ (3) 

The cases where a = 0, 1 are special and the corresponding if-norm (the 
superscript and subscript '0' are omitted in this case) and if^-norm are known 
in the literature as the energy and enstrophy norm, respectively. A geometric 
constraint, referred to as the Poincare inequality, is 



Ml^p>KMl (4) 



for non-negative a, where Ai = A^j, with ki being the smallest wavenumber, 
is the first (smallest) eigenvalue of A in if. The inequality reverses direction 
for non-positive a. 

The bilinear operator -B(-, •) satisfies 

{Au, B{v, v)) + {Av, B{v, u)) + {Av, B{u, v)) = (5) 

for u,v E H^, and 

{u,B{v,w)) = -{w,B{v,u)) (6) 

for u,w & and v E H. These identities arise by virtue of the non-divergent 
and periodic properties of the velocity field. In particular, we have 

{u,B{u,u))^0, (7) 
(Au,B(u,u))^0. (8) 

We shall collectively refer to the above identities as the orthogonality proper- 
ties of the nonlinear term. In the absence of forcing and dissipation, they give 
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rise to conservation of energy and enstrophy. These conservation laws form 
the foundation for the idea of the inverse cascade of energy and the direct cas- 
cade of enstrophy, and the associated concepts of inertial ranges, as mentioned 
earher and reviewed in the next section. 

A well-known mathematical fact is the existence of a bounded finite-dimensional 
global attractor for the traditional 2D Navier-Stokes equations in a finite do- 
main (as considered here). The research performed on this subject constitutes 
a rich literature. The book of Temam [26] provides a good treatment of the 
subject and a complete list of references. Some articles suitable for quick ref- 
erence are Constantin et al. [7,8] and Ziane [27]. We shall not attempt to 
demonstrate the existence of such an attractor for the 2D Navier-Stokes sys- 
tem (driven by a time-independent forcing) with various degrees of viscosity 
(or combinations of them, as presently considered) but take its existence to 
be given. 

In this study / is assumed to be monoscale-like, as defined in the Introduc- 
tion. This means that the ratio of enstrophy to energy injection, which we 
denote by A, is bounded a priori within a certain range and in particular from 
above, for all u. The time-independent monoscale body force / G H{Xs), for 
an eigenvalue As, considered by Constantin et al. [6] is a special case with 
A = As, although it is a rather peculiar special case, for two reasons. First, 
the instantaneous energy and enstrophy injection are not constrained to be 
non-negativeQ, although they clearly must be so in the time mean since the 
dissipation of both quantities is non-negative. This means that the condition 

^ For this reason, a monoscale body force is not, strictly speaking, a special case of 
monoscale-like forcing, even though our results apply to both cases. 
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A = As does not necessarily extend, even approximately, to the case of a 
constant external forcing applied over a range of wavenumbers; such a forc- 
ing is not "monoscale-like" in the sense we require (and as is often assumed 
in the KLB theory). Second, there exists an exact stationary solution (given 
by (26)). Because of this, there appears to be a belief that the flow in this 
case collapses onto the stationary solution and is not turbulent. However, for 
sufficiently strong forcing the forced scales will be unstable to nonlinear inter- 
actions, and numerical simulations (J.C. Bowman, personal communication, 
2001) indeed confirm that there is nothing pathological about the case of a 
constant monoscale forcing: it develops a full spectrum, albeit one constrained 
by the inequalities discussed here. 

The general case of a monoscale-like forcing applied over a finite range of 
wavenumbers, as often imagined in the KLB theory, can be realized in several 
ways. One example (Shepherd [22]) is f — Z^^- cP(Afc)ii/ ||P(Ajk)ii||^ where 
c > 0, P(Afc) is the projection onto H{Xk), \\-\\ is the energy norm, and the 
summation is over the restricted set of wavenumbers K. For this forcing the 
energy and enstrophy injection are both constant in time and A is the mean 
of Xk over K. Another example is / = J2k cP{Xk)u. Although the energy and 
enstrophy injection are now variable in time, both are positive definite and A 
lies within the range [Amin, Amax], where Amin and Amax are, respectively, the 
minimum and maximum Xk in K. This commonly used forcing is known as 
"instability" forcing, for obvious reasons; in the geophysical context it is used 
in 2D turbulence to mimic forcing of the barotropic component of the fiow by 
baroclinic instability (Lilly [16]; Basdevant et al. [1]). In both these examples, 
/ depends on u. Another commonly used forcing is white-noise forcing over 
K (Lilly [16]). This is harder to control a priori, but in practice gives positive 
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enstrophy and energy injection with A generally fluctuating within the range 
[Amin, Amax]- For this case, the condition of monoscale-like forcing would need 
to be verified a posteriori, and would apply at best only in a time-mean sense. 
The results derived here would then apply only in the time mean. 

We now introduce some terminology and derive some preliminary estimates 
which are employed later. In the theory of turbulence the characteristic wavenum- 
ber k defined hy k — \\u\\l^2 / ll'^ll^ often studied. However, we find it more 
natural and convenient in the context of this work to define the parameter 

A=M, (9) 

\\u\\ 

the square root of which has the dimension of a wavenumber. Moreover, let 
Ae/j, and A^^ be defined (for 7^ 0) by 

Ai..= (jj^j , (10) 

/|| ||2 \ V/i 

"^^"^(iRf) ' ^^^^ 

so that 2i/^A^^ and 2i/^A^^ are, respectively, the (instantaneous) energy dissi- 
pation and enstrophy dissipation rates due to Vf^A^^. The values of A^^^ and A^^ 
therefore indicate where, in wavenumber space, viscosity primarily operates in 
the dissipation of energy and enstrophy, respectively, for the dissipation oper- 
ator of degree /i. Note that A^i = A; we will use the latter symbol exclusively 
in what follows. 

We now derive several fundamental inequalities relating these parameters. In 
the following, a, P, 7, and A are real numbers, and (f) is an appropriate func- 
tion so that all norms involved are well-defined. First, we have the following 
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interpolation- type inequality which can be shown by Holder's inequality: 



< 







a+7 117^117 



(12) 



for a > /3 > 0. Note that the inequahty reverses direction for a < /3 < 0. 
Second, let us define 



G{a, (3, X, 0) - - Af) Wn^k 



(13) 



for A > 0. It is obvious from (13) that G(q;,/3,A, 0) is positive (negative) if 
and only if a/? > (a/? < 0) and <p ^ H{X) when A is an eigenvalue of A. This 
implies that when ^ 0, G{a, /3, A, 0) = if and only if A is an eigenvalue 
of A and e H{X). This important feature of G{a, P, X,(f)) is used in some 
arguments of section 4. By rearranging terms we obtain 



A" 



2 >/3 



(14) 



It is then easy to show that 



G\a,i3, 



2/a 



2 II ^||2 

a+/3 
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1 2/a' 



from which it follows that 



(15) 



a+l3 



> 



a IIV^II/3 for "/^ ^ 0' 

p for < 0. 



(16) 
(17) 



The ratio of Ag^ to A^^/ is given by 



A 



En 



^En' 



I i|i/^ II ul/^l'-l/^l' 
\u\\, \\u\\ 

I II II 



\u 



(1^ 



which, by applying (12) with a = /i, (3 = /i' , and 7 = implies 
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Ae/.>Aem'' for// > //' > 0, (19) 
Aem<Aem'> for a* < a*' < 0. (20) 

By applying (12) to the ratio of A^^/A^^/ we also obtain the same inequalities 
for and h-z^'- The ratio of A^^ to A^^ is given by 



, , G(/x,1,A,iz)||m 



|2 



A' 



1+ ^ ; " ■ (21) 



Since G(/i, 1, A, ti) takes the sign of ji (see above), the above ratio indicates the 
simple fact that a viscosity dissipates more enstrophy than energy, while an 
inverse viscosity dissipates more energy than enstrophy Note, however, that 

Azm>A£^, for all At ^ 0. (22) 
We will compare A, A^^ and A^^ with A in subsequent sections. 



3 The KLB theory 



We now review the arguments leading to the KLB theory. In the absence of 
forcing and dissipation, the 2D Navier-Stokes system conserves energy and 
enstrophy due to (7) and (8). There is a consensus that an initial distribution 
of energy cascades towards larger scales, with a downscale cascade of energy 
practically forbidden. Arguments against the downscale cascade of energy were 
advanced by Taylor [24] and Lee [13]. A celebrated "proof of the upscale 
cascade is due to Fj0rtoft [9]. In his argument, Fj0rtoft considered the change 
in energy for three different scales coupled nonlinearly. Because of conservation 
of enstrophy, energy must flow from the intermediate scale to the smaller and 
larger scales, or vice versa. For energy initially on the intermediate scale, it 
was argued that the larger scale acquires most of the transferred energy. In 
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the present setting Fj0rtoft's argument goes as follows. Let I < k < m he the 
three wavenumbers (corresponding to three scales < k^^ < l^^) that are 
involved in the energy transfer. Furthermore, let AE{-) denote the change of 
energy in each scale. It is easy to see from the two conservation laws that 

AE{1) + AE{k) + AE{m) = 0, 
XiAE{l) + XkAE{k) + A^AE(m) =0. 

Hence, 

Am — Xi 

AE{m) = - ^^~\' AE{k). 

Xm " A/ 

It is argued in [9] that if one takes m = 2k = 4/, for example, then (Am = 
4Afc = 16A0 

AE{m) 

Therefore, changes in the kinetic energy are distributed in the ratio 4:1 on the 
components with the double and half scale, respectively, if no other compo- 
nents are involved in the energy transformation. 

There are a number of problems with this argument. First, the inviscid system 
is time reversible, so an upscale cascade of energy cannot be established by 
the above argument alone. Second, the three scales involved in the energy 
transformation must satisfy the triad requirement |^ and not all choices of 
^ Although the geometry considered in [9] is a sphere, interactions in the form of 
triads are to be observed [23]. For an interacting triad of similar scales where one 
wishes to have the smallest and largest scales as far away from the intermediate 
scale as possible, the scale ratio of 1:2:3 is a better approximation than 1:2:4. For 
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the interacting triads give AE{1) / AE{in) > 1 [19]. Nevertheless, Merilees and 
Warn [19] find that an initial spectral peak spreads out with most of the 
energy going upscale — a result very well confirmed numerically. In particular, 
they find that given an intermediate wavenumber k, roughly 70% of interacting 
triads exchange more energy with lower wavenumbers, while roughly 60% of 
interacting triads exchange more enstrophy with higher wavenumbers. Given 
a monoscale forcing at a wavenumber s, then, it is tempting to arrive at the 
picture of the dual cascade of energy to larger scales, and of enstrophy to 
smaller scales. However, in forced-dissipative turbulence with a full spectrum, 
the validity of this picture is not at all self-evident. 

In particular, the KLB theory envisages an enstrophy dissipation range for 
k > kj, with kjj ^ A^/^ (high Reynolds number), so that the forcing and dissi- 
pation scales are well separated. If molecular viscosity is the only dissipation 
mechanism, then the dissipation of enstrophy is at least A;^ times the dissipa- 
tion of energy. On the other hand, the forcing of enstrophy is only A times 
the forcing of energy. This would suggest that forced-dissipative equilibrium 
is unrealizable for high- Reynolds- number 2D Navier-Stokes turbulence. Yet, 
Constantin et al. [7,8] have proven the existence of a global attractor for this 
system for the case of time-independent forcing. Unless that case is patholog- 
ical, it follows that the assumption of an enstrophy dissipation range confined 
to A; > A;,^ ^ A^/^ must be wrong. In order to deduce the directions of energy 
and enstrophy cascades, one needs to know the spectral distribution of energy 
and enstrophy dissipation, and these are not preordained. 



such a case, the ratio of the upscale cascade energy to the downscale cascade energy 
is approximately 5:3. 
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4 Asymptotic behaviour: Dynamical constraint 



We first derive the governing equation for the evolution of A in the absence of 
a forcing term. Taking the scalar product of (1) with u and Au, respectively, 
we obtain the governing equations for the decay of the energy and enstrophy 
for / = (note that the nonlinear term vanishes in both cases due to the 
orthogonality properties) : 



"^f + 2E^mII< = 0' (23) 



dt 

^^\\u\\l + 2j:uju\\l^^ = 0. (24) 



Note that the subscript n has been dropped from the sum as there is 
no risk of confusion. Taking the time derivative of the expression for A and 
substituting (23) and (24) we obtain 



2 II ||2 



dA_ ^ Fill fIU - I|m|| \\u 



dt 



\u\ 



|4 



= -2 2:.„%^. (25) 



\u\ 



A couple of remarks are in order: 

Remark. If Ekman drag (/x = 0) is the only dissipation mechanism involved, 
then dA/dt — because G(0, 1, A, u) — 0. In this special case both enstrophy 
and energy are dissipated at the same rate 2i/o- Since A remains constant, the 
distribution of energy and enstrophy should then be dramatically different 
from the viscous (// > 0) cases. 

Remark. A viscosity causes A to decrease monotonically since G(/i, 1, A,u) > 
for all /X > 0, while an inverse viscosity has the opposite effect. If a single 
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viscosity mechanism is involved, the rate of decay of A is greater for a more 
spread out spectrum as compared with a sharper spectrum having the same 
energy and enstrophy. This can be seen by the fact that l,A,u) is greater 
in the former case. 

We next consider the time- independent monoscale forcing / e H{Xg), where 
As > Ai, considered by Constantin et al. [6]. For the arguments to follow we 
note that (1) possesses a stationary solution given by 

^-(E^.KY'f- (26) 

This stationary solution is referred to as the primary stationary solution. A 
stationary solution other than u (if it exists) is identified as a secondary sta- 
tionary solution. The existence of such a solution for the traditional 2D Navier- 
Stokes system under suitable conditions is demonstrated in [10]. The energy 
and enstrophy evolve according to 

^^rf + E^.lK = K/)> (27) 
^^IHI? + E^mIHI?+, = (^^,/)- (28) 

Since {Au, f) — Xs{u, f) in this case we can multiply the energy equation (27) 
by Xg and subtract it from the enstrophy equation (28) to obtain 

{hWl - ll^ll'} + E^M - \Ml}^0. (29) 

This equation can be rewritten in terms of the function G{ii, 1, Xs, u) as 



Ift+T. ^^Ki = - E 1' ^s. u), (30) 

1 1 2 1 1 2 

where ^ = \\u\\-^^ — Xs \\u\\ . Now we require that jj, be non-negative, i.e. no 
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inverse viscosity is allowed, so that the right-hand side of (30) is non-positive. 
This means that any positive value of ^ will be dissipated away, while a non- 
positive value will evolve but remain non-positive for all time. To see this 
exphcitly, consider the formal solution of (30), for all t > to, given below: 

e(t)=exp{-2X:i^MA^4 X (31) 

where C,o = ^{t = to)- It is easy to see that if ^ then ^{t) < for all time. 
On the other hand, if > then ^{t) becomes non-positive for all i > T 
where T is the solution of 

T 

Co exp {2 ^ i/^A^io} = 2 y exp {2 5^ i/^A^i'} ^ u^G(ii, 1, A„ u)dt'. (32) 

to 

Equation (32) always has a finite solution for T unless ^' -^s' ^) 

creases exponentially in time more rapidly than exp{—2J2^nXs^}- Since the 
function G{ii, 1, A^, u) vanishes if and only if u e H{Xs), this can occur only if u 
asymptotically approaches (with respect to the various norms) the eigenspace 
H{Xs). Now H{Xs) is part of the stable manifold oi u — (Z^ i^^uA^)"^/; this 
can be seen from the fact that B{u,u) — for all u e H{Xk), \/k, so H{Xs) 
is invariant and every trajectory on it exponentially converges to u at the 
rate 2 J2 ^n^t- K follows that any trajectories which asymptotically approach 
H{Xs) also converge to u. Therefore, the cases for which a finite solution for T 
is questionable turn out to be restricted to trajectories on the stable manifold 
of u. 

The argument in the last paragraph indicates that for a general trajectory 
not on the stable manifold of -u, ^ acquires a non-positive value in a finite 
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time. Within the stable manifold of u, only u lies on the global attractor and 
^ = for M = M. But before drawing a final conclusion for C, on the global 
attractor it is necessary to secure the non-positiveness of ^ on homoclinic 
trajectories of u which, if they exist, are part of the global attractor. Recall 
that a homoclinic trajectory (also known as a homoclinic orbit) of a stationary 
solution is one that emanates from the stationary solution and terminates on 
that same stationary solution. It is a special limit cycleQ with infinite period. 
It can be seen from (31) that ^ does not acquire a positive value on any 
homoclinic orbit of u during its infinitely long journey, if such an orbit exists. 
It would, rather, start from u with = 0, evolve with some negative value for 
^, and then ^ would increase to terminate at u with = again. 

Therefore, it can be concluded that on the global attractor of the Navier- 
Stokes system, the energy and enstrophy satisfy 

< As . (33) 

This inequality will be referred to as the dynamical constraint. It is noted, for 
the sake of completeness, that the argument leading to (33) also implies that u 
is the only point on the global attractor where the equality occurs. Inequality 
(33), together with the Poincare inequality, gives 

Ai < < As . (34) 

^ Note that a regular limit cycle which would intersect H{Xs) does not exist because 
H(Xs) is part of the stable manifold of u. Alternatively, the existence of such a cycle 
would imply from (31) that Jq° exp{2 ^ u^Xgt'} ^ v^G{^, 1, A^, u)dt' = 0, where Tq 
is the period of the cycle. But this would require the cycle to be entirely in H{Xs), 
which would further reduce the cycle to u. 
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An immediate corollary of the dynamical constraint is a similar constraint 
on the enstrophy of a secondary stationary solution if one exists. In fact, 
(33) becomes a strict inequality for any secondary stationary solution. More 
interpretation of (33) is given in the next section. 

Remark. The addition of an inverse viscosity to a system which possesses 
a global attractor does not jeopardize its existence. However, the spectral 
distribution of energy and enstrophy on the attractor does not necessarily 
obey the dynamical constraint in this case. 

Remark. The requirement of a time-independent / is merely for the sake of 
securing the existence of a global attractor. Equation (33) is a constraint on 
the long-time behaviour of the dynamics whether or not the monoscale forcing 
/ G H{Xs) is time-independent. 

We now extend the above result to the case of a more general monoscale-like 
forcing as defined in the Introduction. It is easy to see that when A is constant 
the above analysis applies with A in place of A^. (In general, there will not 
exist a stationary solution as in the case of a monoscale forcing.) For variable 
A, we have the a priori inequalities (by hypothesis) 



where Amin and Amax are, respectively, the minimum and maximum A^ of the 
restricted set K of wavenumbers through which energy and enstrophy are 
injected. Note that with variable A, we find it necessary to require the condi- 
tion (uj) > at all times. Then (27), (28), and (35) lead to (30) with the 
equality replaced by < and with A^ replaced by Amax- It follows that ^{t) is 
bounded from above by zero, for sufficiently large time, and thus the dynam- 





(35) 
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ical constraint (33) must hold on the attractor (if it exists) with replaced 

by Amax- 

Finally, if (35) only holds in a time-mean sense — note that the time-mean 
energy injection {{u, f)) > necessarily since the energy dissipation is non- 
negative — then we may take the time mean of (30) (with A^ replaced by Amax 
and the equality replaced by <) to obtain {^) < 0, equivalently 



and the dynamical constraint holds in the time mean. 
5 Dissipation rates and spectral distributions 

In this section we examine the time-mean dissipation rates of energy and en- 
strophy and possible spectral slopes of energy for three combinations of gen- 
erahzed viscosity: (i) molecular viscosity alone, (ii) molecular viscosity and 
Ekman drag, and (iii) molecular viscosity and inverse viscosity. Since we are 
assessing the vahdity of the KLB theory, which applies to statistically station- 
ary turbulence, we may assume that all time-mean quantities exist. 

For simplicity we first treat the case of monoscale forcing as in the last section, 
leaving the question of the applicability of the results to the more general 
monoscale-hke forcing until the end of the section. For the first two cases, 
the results from the last section indicate that both energy and its dissipation 
are confined to scales no smaller than the injection scale. In particular, the 
dynamical constraint gives 




(36) 
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A = 




< A, 



(37) 



This means that in both cases the spectrum adjusts so that the dissipation of 
energy takes place on scales no smaller than the forcing scale. Moreover, the 
energy dissipation rates are bounded from above by 2i>iXs for case (i) and by 
21/1A5 + 21/0 for case (ii). 

We know that viscosity acts on smaller scales for enstrophy than it does for 
energy. The question is where, in wavenumber space, the enstrophy dissipation 
takes place. To answer this question we resort to (29). Since both energy and 
enstrophy are bounded from above on the global attractor, taking the time 
mean of (29) leads to 



This equation is referred to as the balance equation. Attempts to constrain 
the time-mean enstrophy dissipation rate and the scaling law of the energy 
spectrum constitute the remainder of this section. 

Case (i): The balance equation reads in this case 



This well-known result corresponds to (6) of Constantin et al. [6]. It implies 
that the time-mean enstrophy dissipation rate is given by 2z/iAs, and that the 




(38) 



{Ml) = K{\\u\\l), 



(39) 



so 




(40) 



22 



enstrophy dissipation is concentrated around the forcing scale. It follows that 
there can be no direct enstrophy cascade. Moreover, this constraint gives rise 
to a very steep spectral slope for k > s. Let E{k) be defined by 



E{k)^{\\P{X,)uf) 



so that the total energy £ is given by 



oo 



£^ / E{k)dk, 







where the discrete wavenumber has been changed to a continuous one for the 
sake of convenience. Now suppose that E{k) oc k~^ over a range [kg, ki] where 
s < ko <^ ki; the dissipation of enstrophy over that range, Dz, is then given 



Equation (40) requires that Dz be dominated by the large scales. As a result, 
the satisfactory values of 6 are confined to S > 5. {S — 5 is not acceptable.) 

Remark. A scaling law of the energy spectrum k"^ with 5 > 5 is much steeper 
than the classical inertial range power law (see Kraichnan [11,12]). Hence, 
the result of this section suggests that the latter scaling law is unreahzable 
in 2D Navier-Stokes turbulence for a constant monoscale forcing. This result 
was derived by Constantin et al. [6]. 

Remark. For a viscosity of degree /i > instead of molecular viscosity, which 
includes Ekman drag (// = 0), the satisfactory values of 5 are confined to 
5 > 3 + 2/x. 



by 




(41) 
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Remark. Although an inverse energy cascade cannot be excluded, we can rule 
out the classical k~^^^ power law [11] for any such range by a similar argument. 
Suppose that E{k) oc k~^ over a range [kg, ki] where ko ^ ki < s; the 
dissipation of energy over that range, De, is then given by 

ki 

DecxJ k'-^dk = ^ {kl-^ - kl-") [7 3]. (42) 

ko 

An inverse cascade of energy requires that De he dominated by the small-/c 
end of the range, which requires 7 > 3. (7 = 3 is not acceptable.) Therefore, 
the k~^^^ power law for 2D Navier-Stokes turbulence is incompatible with 
an inverse energy cascading range. This result is clearly independent of the 
nature of the forcing. 



Case (ii): The balance equation reads in this case 

{{\\u\\l) - Xs{\\u\\l)} + 1.0 {{\\u\\l) - Xs{\\uf)} = 0, (43) 



1^1 



so 



uiAzi + 1^0 = uiXs + > i^iXs + uq. (44) 

This implies that the time-mean enstrophy dissipation rate is bounded from 
below by 2i/iAs-|-2i/o- Equivalently, we have Azi > A^. In other words, Ekman 
drag allows the spectrum to adjust toward the small scales (as compared with 
case (i) for which A^i = A^), so that the enstrophy dissipation can occur at 
scales smaller than the forcing scale. To see the extent of such an adjustment 
we rewrite the equality in (44) as 
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The ratio A-z\/\ > 1 because A^/A > 1 (strictly speaking we need to exclude 
the primary stationary solution rt, the only point on the global attractor for 
which As/A = 1). The question is whether Kzi/\s can be much greater than 
unity for a system with vq ^ i^iA^, where the large-scale dissipation is dom- 
inated by Ekman drag. This is not obvious since we expect that A^/A 1 
as i^q/{i^i\s) ~^ because the dynamical constraint is A = A^ when Ekman 
drag alone is responsible for dissipation. However, we do not know exactly 
how this limit is approached as fo/ifiXg) ~^ oo- Nevertheless, it appears that 
in principle a direct enstrophy cascade could be achieved in this system. 

We now examine implications for scaling of the energy spectrum. For an energy 
cascading range /cq <C /ci < s, the energy dissipation is given, with the same 
notation as before, by 

ki 

De(xJ {iyok~'^ + yik^"^) dA; 

feo 

= -^(k\---kl-'^) + ^(kl-'^-kl-'^) [7v^l,3]. (46) 

We have already assumed that z/q S> i^iXg so that Ekman drag dominates on 
the large scales, and therefore the energy dissipation occurs primarily on the 
large scales provided 7 > 1. On the other hand, we want negligible dissipation 
of enstrophy on the large scales, and this requires 7 < 3. Thus 1 < 7 < 3, 
which includes the 7 = 5/3 power law of the KLB theory. 

For an enstrophy-cascading range s < <^ ki, the enstrophy dissipation is 
given by 
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fcl 

Dzoij (uok^-^ + uik^-^) dk 

ko 

= 33l(^i''-^o-V^(^i"'-^o"') [^^3,5]. (47) 

It is seen that enstrophy dissipation will occur primarily on the small scales 
provided 5 < 5. {5 — 5 is not acceptable.) On the other hand, we want 
negligible dissipation of energy on the small scales, and this requires 5 > 3. 
(5 = 3 is not acceptable.) Thus 3 < 5 < 5, and the 5 = 3 power law of 
the KLB theory (with or without the logarithmic correction of [12]) is not 
permitted. In fact, it is in any case excluded (as in case (i)) by the dynamical 
constraint (37). 



Case (iii): Similar to case (ii) we have 
where 



A^,^^, A^,^^^^. (49) 



Unlike case (ii) where only one ratio on the right-hand side of (45), namely 
A5/A, depends on the actual spectrum, all three ratios on the right-hand side of 
(48) are spectrum dependent. Among them we expect A^^/A^;^ to be smaller 
than unity for negative /i (see section 2). Now the ratio A^/A is not always 
greater than unity, and depends on the outcome of the competition between 
molecular viscosity and inverse viscosity. Nevertheless, for suitable values of 
and /I a dual cascade certainly appears to be possible. In this case, energy 
cascades upscale and is dissipated by inverse viscosity, while enstrophy cas- 
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cades downscale and is dissipated by molecular viscosity. Of course, inverse 
viscosity has no physical basis and is employed purely for numerical reasons. 

We now extend the results of this section to the more general case of a 
monoscale-like forcing. It is only necessary to do so for cases (i) and (ii) since 
those are where we have derived restrictions. It is easy to see that every- 
thing goes through when A is constant. For variable A, or if the condition of 
monoscale-like forcing holds only in the time mean, the balance equation (38) 
may be replaced by 



The arguments concerning admissible spectral slopes go through without change 
since they do not refer at all to the forcing mechanism; in the definitions of 
the inertial subranges the forcing wavenumber s is replaced with either A^„ 
or Xll^, as appropriate. This leaves only the prohibition of a direct enstrophy 
cascade in case (i) . There one obtains Azi < Amax in place of (40) , which still 
precludes a direct enstrophy cascade. 

6 Concluding remarks 

In this paper we have derived various constraints on the asymptotic be- 
haviour of the 2D Navier-Stokes equations in a finite domain driven by a 
monoscale-like forcing and damped by a general class of dissipation operators. 
By monoscale-like forcing we mean that the ratio of enstrophy to energy in- 
jection lies within the range of the square of the forcing wavenumbers, for all 
possible velocity fields u, which is a common (though not exclusive) scenario 
in the KLB theory (e.g. Kraichnan [11], p. 1421b; Pouquet et al. [21], p. 314; 




(50) 
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and Lesieur [15], p. 291). Note that this crucial property is not guaranteed for 
a constant external forcing, except in the special case of a monoscale forcing. 
The results obtained include constraints on the time-mean enstrophy dissipa- 
tion range, an upper bound on the system's enstrophy in terms of the energy 
and the forcing scale, and bounds on possible scaling laws for the energy spec- 
trum. If the condition of monoscale-like forcing holds only in the time mean, 
then the upper bound on the enstrophy holds in the time mean rather than 
asymptotically. The validity of the dual cascade picture is examined by de- 
termining where, in wavenumber space, dissipation of energy and enstrophy 
occurs. 

For 2D Navier-Stokes turbulence, with molecular viscosity the only form of 
dissipation, no direct enstrophy cascade is permitted; rather, enstrophy dis- 
sipation is required to occur in the vicinity of the forcing scale. Although a 
reverse energy cascade is permitted (not to say that it occurs), it is shown 
to be incompatible with the —5/3 power-law scaling of KLB theory for the 
energy spectrum in such a range (as noted previously by Constantin et al. [6]). 

Ekman drag has often been employed in numerical simulations of 2D turbu- 
lence to allow a second dissipation channel. Furthermore, Ekman drag is a 
reasonable representation of a frictional planetary boundary layer in the geo- 
physical context. We show rigorously how the introduction of Ekman drag to- 
gether with viscosity breaks the strong constraints on the spectral distribution 
of enstrophy dissipation, and thereby allows the possibility of a dual cascade. 
However, the rigorous bound on the enstrophy precludes the existence of the 
—3 power-law scaling of KLB theory (with or without a logarithmic correc- 
tion) for the energy spectrum in an enstrophy-cascading range. Instead, the 
power law is shown to be between —3 and —5. This is consistent with the 
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published results of numerical simulations. 

It is shown that the use of an inverse viscosity together with viscosity allows 
the possibility of KLB power-law scaling in the enstrophy-cascading range. 
Indeed, the numerical simulations of Borue [4], and the more recent high- 
resolution numerical simulations of Lindborg and Alvelius [17], which both 
claim to exhibit this scaling, employ an inverse viscosity. Our results suggest 
that this is not a coincidence. 

It should be emphasized that our derived constraints on power laws for the en- 
ergy spectrum apply for any spectrally localized forcing, not just for monoscale- 
like forcing. 

The results derived in this paper are believed to be particular to 2D turbulence, 
because they all rely on the conservation of both energy and enstrophy in 
nonlinear interactions, as expressed through (27) and (28) and subsequent 
key relations such as (29), (38), and (50). There is no analogue of (28) in 3D 
turbulence. 

Constantin et al. [6] have shown that for the special case of a spectrally- 
localized constant external forcing, the realization of a direct enstrophy cas- 
cade in forced 2D Navier-Stokes turbulence in a finite domain requires a com- 
bination of energy input and energy removal by the forcing. The generalization 
of this result to our case is that a direct enstrophy cascade within forced 2D 
Navier-Stokes turbulence requires the abandonment of the popular concept of 
monoscale-like forcing (as defined here). Of course, the —5/3 power- law scal- 
ing remains incompatible with an upscale energy cascade, irrespective of the 
forcing mechanism. 
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A general result of this study is that the choice of forcing mechanisms and 
dissipation operators has implications for the spectral distribution of energy 
and enstrophy dissipation, and thus for the possible existence of energy and 
enstrophy cascades. Furthermore, the choice of dissipation operators has im- 
plications for permissible scalings of the energy spectrum. In choosing forcing 
mechanisms and dissipation operators for numerical reasons, one should be 
mindful of these constraints. 
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